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The integral relaxation time describing the thermoactivated escape of a uniaxial quantum spin sys- 
tem interacting with a boson bath is calculated analytically in the whole temperature range. For 
temperatures T much less than the barrier height A?7, the level quantization near the top of the 
barrier and the strong frequency dependence of the one-boson transition probability can lead to the 
regularly spaced deep minima of the thermoactivation rate as a function of the magnetic field applied 
along the z axis. [S1063-651X(97) 13203-2] 
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The problem of the escape rate of a uniaxial magnetic 
particle has remained in the focus of attention since the 
work of Neel [Q , who stressed the role of thermal agita- 
tion. During the last years the interest in this problem 
has increased in view of possible applications to the in- 
formation storage and in connection with the magnetic 
quantum tunneling (MQT; see e.g., Ref. ||^). 

For classically large particles, the thermoactivation es- 
cape rate was first calculated by Brown who derived 
the Fokker-Planck equation for an assembly of parti- 
cles and solved it perturbatively in the high-temperature 
case, Af/ <C T, and with the use of the Kramers 
transition-state method Q for T ^ AC/. In these both 
limiting cases the time dependence of the average mag- 
netization is a single exponential, and the relaxation rate 
of ferromagnetic particles is given by the lowest eigen- 
value Ai of the Fokker-Planck equation. For T ~ AC/ 
the latter is no longer the case, and the best measure 
of the relaxation rate is the integral relaxation time Tint 
determined as the area under the magnetization relax- 
ation curve after a sudden infinitesimal change of the 
longitudinal magnetic field The quantity Tint 

can be calculated analytically in the whole range of tem- 
peratures, and coincides with Ai in the asymptotic 
regions. 

With the miniaturization of the magnetic particle both 
the thermoactivation and the MQT escape rates increase; 
the latter becomes dominant below the crossover tem- 
perature To determined by the interactions noncommut- 
ing with the operator Sz and hence causing the MQT. 
For information storage applications the most important 
are systems with small tunneling interactions and cor- 
respondingly low Tq. In this case the tunneling level 
splittings can be calculated perturbatively P| for the 



arbitrary spin values S, which can be advantageous in 
comparison to the semiclassical instanton method (see, 
e.g., Ref. Q) for nanoscale systems with moderately high 
spin values as the recently synthesized Mn clusters hav- 
ing 5 = 10 and 12 in the ground state 11|. For such 
systems with not too large S and low To the thermoacti- 
vation escape rate can be dominant down to the tempera- 
tures where it changes due to the spin level quantization. 

The purpose of this paper is to calculate the thermoac- 
tivation escape rate of a quantum spin system in terms 
of the integral relaxation time starting from a spin-bath 
Hamiltonian of the type 



n = -HSz - DSl + Tib 
-^l-q(r/S)(at + 
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where the arbitrary vector t/ is for simplicity set to be 
rjx = rjy — rjz = I. If the bath excitations described by 
the operators and a are phonons, then the coupling 
to the bath that is linear in spin variables is prohibited 
by the time-reversal symmetry. This means that modu- 
lations of the crystal field by phonons do not produce a 
fieldlike perturbation on a spin system, and hence r] = 0. 
Thus, in this case it would be better to write quadratic 
terms of the type SaS/^ instead of {rjS) in (|l|) (see, e.g., 
Ref. ||l^). Nevertheless, we will use (|l|) with Tit and 
the couplings V characteristic for phonons since it is the 
most suitable for the first presentation of the method and 
describes the main qualitative features of the relaxation 
of a spin system. Moreover, Hamiltonian (|l|) means the 
direct quantum generalization of the Langevin-/ie/d for- 
malism used by Brown |^ (which is subject to the same 
criticism), and thus provides a link to the known results 
for classical ferromagnetic particles. 

If the spin-bath coupling is weak, the equation of mo- 
tion for the density matrix of the spin system can be 
obtained in the second order of perturbation theory, and 
the diagonal part of it is the well-known system of the ki- 
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netic balance equations for the occupation numbers Nm 
of the spin states \m) 

+ lm-liWm,m-lN,n-l " W,n-l,mNm)- (2) 



It is convenient to introduce the reduced variables 



Here Im = \/S{S + 1) — m{m + 1) are the matrix ele- 
ments of the operators S± , and the " spin-free" transition 
probabilities W are given by 



where 



'^m+l,m = £m+l — £m — —H — D{2m + 1) 



(4) 



are the transition frequencies and £,„ = —Hm — Dm? are 
the spin energy levels, 

q 

+ TlqTr6{ujq — U})] (5) 

is the contribution of the onc-phonon emission and ab- 
sorption processes, and 

W^(2)(^) = \Vpqfnp{ng + l)TrS{Ljp - cj, - cj) (6) 



is that of the two-phonon (Raman scattering) ones. One 
can check that the transition probabilities satisfy the de- 
tailed balance condition: W{uj) = exp(— w/T), 
which ensures the static distribution of the form 

1 ^ 
iV^o) = - exp(-e™/T), Z = exp(-e™/T). (7) 

m=-S 

The quantities W(w) with w > describing transitions 
to upper energy levels become exponentially small at low 
temperatures. For spin-phonon couplings of the type 



(ujpUJq) 



1/2 



(8) 



where il = Mc , M is the unit cell mass, ujq — cq, and 
c is the phonon velocity, the estimation of W{uj) with 
w < with the help of (|) and (|) yields (cf. Ref. p|) 



and 



g2|,.,|3 
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-("It^l + 1) = 



(9) 



020|,TVe8, \lu\ ^Od^T 
^ I elT^/Q\ \u\ « T « Od (10) 

Here 6'_d ^"^gmax is the Debye temperature, = 
r26'|3 ~ h? p^C', and p is the density of the lattice. 
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which are equivalent to those used for the description 
of classical single-domain magnetic particles ^,|[ and 
should be kept constant ii S ^ oo. In this limit the 
transition frequencies uJm+i,m of Eq- (jj) tend to zero 
and, accordingly, the frequency-independent two-phonon 
transition probabilities W^'^^ given by (^0|) govern the re- 
laxation. Since the occupation numbers of the neighbor- 
ing levels in (^) become close to each other, Eq. 
over to the classical Fokker-Planck equation 
extreme quantum case is realized for a three-level system 
with a barrier {S = 1), for which the relaxation between 
the two lowest levels through the highest one (the so- 
called "resonance fluorescence" ) has at low temperatures 
an exponentially small rate F ~ exp(— A/T) p3[ . 

The relaxation of any initial state described by the 
system of the first-order differential equations (^ is de- 
scribed in general, and particularly at intermediate tem- 
peratures, by 25* exponentials of the type exp(— A^t), 
where are the corresponding eigenvalues. In such sit- 
uations the best measure of the relaxation rate is the 
integral relaxation time Tint determined as the area un- 
der the magnetization relaxation curve after a sudden 
infinitesimal change of the applied field H: 



Tint 



m^((X)) - mj;(0)' 



(12) 



Here = {Sz)/S is given at equilibrium by — 
B{^,a) = dlnZ/d^. One can calculate Tint analytically 
at arbitrary temperatures from equations (|^) producing 
the low-frequency expansion of the linear longitudinal dy- 
namic susceptibility 



Xz(w) = Xz(l + iwTint + ...). 



(13) 



Here Xz = drriz/ d{SH) = B' /T is the static suscepti- 
bility, and B' = dB{S^,a)/d£,. Taking into account that 
the alternating field AHz{t) — Ai/^o exp(— iwt) modu- 
lates the transition frequencies (|^), using the detailed 
balance condition and introducing N„i = Nm\i + Qm) 
with Q,n = Qmi'-^)SA.Hz{t)/T one comes in the linear 
approximation in AH^ {t) to the system of equations 

m+1 Qm 

+ iujqm = {li,W,n+l,,n - W^m-l,m)/5'. (14) 

The susceptibility Xz('j-') can now be written as 



Xz(w) 



TS 



J2 mN^kmii^). 



(15) 



The second-order finite-difference equation (|lj) should 
be solved perturbatively in uj; this can be done analyti- 
cally since the first line of (p^) c ontains only the differ- 
ences of Qm, and the order of (|14|) can thus be lowered to 
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1. In the static limit. 

,(0) 



y) reduces to q!;;:'^^ - ql^^ = 1/S 
with the solution qin' = m/ S — B. In first order in uj one 
can introduce = Z^Wm+i,™ exp(-e™/r)[g^^|i-gli^], 
satisfying = -i^g™'' exp(-£m/r). Finding p,„ 
and then qm\ and using jisl ) and (p^, one obtains the 



final result 
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Tint 
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(16) 



where A'^^^ is given by and 



E 

fe=-S 
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N, 



(0) 



(17) 



The formulas above are valid at all temperatures, and 
they are a direct quantum generalization of the classical 
results of Ref . , which are recovered in the limit S 
00. At high temperatures where a, f <^ 1 and VF*^^-* <^ 
W^'^\ the calculation in (p7|) and (O) yields 




(18) 



with An = 2iy(2) (pf^ j|g£ |jn^ At a,^ - 1 (i.e., 
T ~ A[/) the sums in ( |l7| ) and (^ should be performed 
numerically. The recent numerical calculations of Ref. 
for the classical model have shown that, in the unbiased 
case, = H = 0, the difference between and the 
lowest eigenvalue of the Fokker-Planck equation Ai does 
not exceed 1.2% in the whole range of temperatures. Ai 
remained in the focus of interest since the work of Brown 
but for T ^ AU it has no direct physical meaning 
and cannot be represented by a closed analytical formula. 
One can expect that the quantities T[~j and Ai are even 
closer to each other in the unbiased quantum case, since 
the difference between them disappears for S = 1/2. 

At low temperatures a ^ 1, the thermoactivation es- 
cape rate of the particle becomes exponentially small. 
Here, as was recently discovered in Ref. Q], t.~1 Ai 
for sufficiently strong bias {h > 0.2). This effect was phys- 
ically explained ||^ as resulting from the depletion of the 
upper potential well and the competition in the integral 
relaxation time Tint of the overbarrier thermoactivation 
with the rate Ai and the fast relaxation inside the lower 
well. These two relaxation mechanisms can be analyt- 
ically separated for S* 3> 1 since the summand of ( |l^ ) 
consists for a, ^ 3> 1 of two peaks centered at the barrier 
top and in the lower well (m S) . The barrier contribu- 
tion Tint,B can be related to Ai taking into account the 
depletion effect in the small bias case, ft, ^ 1, which 
will be considered henceforth, one obtains 



Ai^ =Tint,BB'(e,a)cosh2e■ 



(19) 



(0) 



For ^ = one has i?' = 1 and Ai = t^^^ . Since iV^."' m 
( |l7| ) is strongly peaked at low temperatures in the wells 



{k iS), and is small elsewhere, the function is 
independent on m and given by $m = 1/(2 cosh^^) in 
the main part of the phase space, including the barrier 
region. Then ( pj|) and (|l6| ) can be combined to 



Ai 



ASjS + 1) cosh^ g 



^ exp(£m/r) 



m— — 00 



m+l,m 



where the exponential factor ^ ^/^k^'^ in (^OD cuts the 



(20) 



sum actually at m Sa^^/'^ < S [see (|) and (|T^)]. The 
partition function Z in (^) for ft ^ 1 is given by 



(5/0)6" cosh SD <t: S^D 



2e" cosh^. 



(21) 



where SD is the level spacing in the wells. 

The sum in ( pO| ) depends on the relation between tem- 
perature T and the level spacing near the top of the bar- 
rier '--^ D, as well as on the role played by the one- and 
two-phonon probabilities W^^^ and W^'^\ For not too 
low temperatures this sum can be approximated by the 
integral over x = m/S, and the quantity W^^'^ of (^) can 
be in this case represented in the form 



(1) 



D2 
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(22) 



whereas W^^'^ oc is given for T <C 9d by the second 
line of (|l^). One can see that the integral over x in 
(pO|) can be cut either by the exponential function at 
Ax ~ a~^/2 or by the denominator W{u!m+i,m) at Ax ~ 
y^W^^^yjW^ . Crossover between these two regimes 
occurs at the temperature 



1/5 



2/5 



(^12 « 0d). (23) 



Taking ( |2l| ) into account, one obtains 

Ai = A{ exp[-a(l -I- hf] + exp[-a(l - hf]}, (24) 
where ft <C 1, and the prefactor A is given by 



21^(2)^-1/2^3/2 



^^514^(2)^-1/2^1/2 ri2«r«57^ 



The first of these expressions coincides with that of 
Brown in the high-barrier limit, if one introduces 
An = 211^(2) TJie temperature dependences of A read 
T"/2, T^, Ti3/2^ and T*, respectively. 

The continuous approximation above is only valid if 
Ax - S-^y^W^^yW^^^ > S-'\ or Am > 1. This con- 
dition sets one more charakteristic temperature 



1/3 



1/3 



(26) 
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below which the level quantization becomes essential. In 
the temperature interval D <C T <Tg the sum in ( |20| ) 
converges at Am ^ 1, whereas the exponential factors 
can still be neglected. In this case Ai shows a sinu- 
soidal dependence on the longitudinal magnetic field in 
the weak- field region H ^ D, the amplitude of which is 
exponentially small for W^'^'^ > W^^"'^^ but becomes great 
for ^ ^(1) i-jjg j.j^jjgg jy ^ T<Tq,SD one 

obtains formula (pj), with the prefactor given by 



^ 2D 2 



wW + vi^(2) I ^ (27) 



This result shows deep minima of the thcrmoactivation 
escape rate Ai for H/{2D) = ±1/2, ±3/2, ... for in- 
teger and for H/{2D) = 0,±1,±2 ... for S* half-integer. 
At such fields two levels near the barrier top become de- 
generate, and the leading contribution to the transition 
probability between them, disappears. 

In the extreme quantum temperature region T <C -D, 
the sum in ( |20| ) is again determined by the exponential 
factor. Here, however, only one or maximally two terms 
of the sum contribute to (^0|). If there is an energy level 
with m = Tninax just at the top of the barrier, like for 
S integer and H — (mmax — 0), then equal contribu- 



tions to ( |2C| ) come from the terms with m = mmax and 
m = TOmax ~ 1 , which describe transitions between three 
upper levels with m — rTimax, 'm-ma.x ± 1. If there are two 
degenerate levels to the right and left from the barrier 
top, like for S half-integer and H = 0, then the leading 
contribution is due to the term in ( |20| ) describing tran- 
sitions between these two levels. The general result is 
given by formula with 



A^S{S + 1) cxp(-Ae„„^jT) 
1 1 



W, 



w, 



(28) 



where m^ax = ~H/{2D) + F, F = T[H/{2D) + S + 
1/2] — 1/2, Aeto^^^^ — DF^ is the mismatch between 
£"imax aiid the barrier top, and J-{X) is the fractional 
part of X. The transition probabilities W in (|2^) are 
dominated by W^^'' oc [see (||)] with WT?x„ax±i,nimax — 
— Z?(l ± 2F), and the escape rate Ai shows qualitatively 
the same magnetic field dependence as that given by (27). 
Note that in this low-temperature limit the prefactor A 
is temperature-independent, excluding the narrow field 
regions where it is determined by VF^^^ and is very small. 

Summarizing, the thcrmoactivation escape rate of a 
quantum ferromagnetic particle was calculated micro- 
scopically in the whole temperature range, allowing for 
the frequency dependence of the transition probabilities 
and for the quantization of the energy levels. Even in 
the low-temperature range, where the escape rate is ex- 
ponentially small, the situation is determined by several 
characteristic energies and there are rather many limiting 
cases for the prefactor A in (|2^), which can be difficult 
to observe if the corresponding temperature intervals are 



not wide enough. Not trying to give numerical estima- 
tions, since the Hamiltonian (|l|) is only a schematic one 
suitable mainly for the presentation of the method and 
for a qualitative analysis, we make only some general re- 
marks about the results obtained. First, the greatest va- 
riety of different temperature intervals can be realized for 
particles containing a macroscopically large number N of 
magnetic ions and having an effective spin S'cff ^ N ^ 1. 
However, only the classical case [the first line of (^] can 
be practically observed, whereas in other cases Ai is un- 
measurably small due to the too large values of a. 

Better candidates for searching for the nonclassical 
thcrmoactivation rates predicted here are nanoscale sys- 
tems such as Mni2 clusters with 5 = 10 having a strong 
uniaxial anisotropy (At/ = S^D = 61 K) 0. These 
clusters show a superparamagnetic behavior, and for 
2 K < T < 8 K the prefactor A in (||) is tempera- 
ture independent. The latter corresponds to the extreme 
quantum case (^), which could, however, be expected 
only for T<D ^ 0.6 K. In the main range of temper- 
atures T < S^D = 61 K, one cannot expect one of the 
pure limiting forms of A, and should resort to using ( [20| ) 
because 5* = 10 is not large enough. Nevertheless, one 
qualitative feature always remains: If the relaxation is 
governed by the two-phonon processes, one can expect 
a strong temperature dependence of A, and, if the one- 
phonon processes are dominant, then the temperature 
dependence of A is weak or absent, but there should be a 
strong dependence on the magnetic field of the type ([2^). 
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